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Suspensions of copper nanowires are emerging as new electronic inks for next-generation flexible electronics. Using a novel 
surface acoustic wave driven extensional flow technique we are able to perform currently lacking analysis of these suspensions 
and their complex buffer. We observe extensional viscosities from 3 mPa-s (1 mPa-s shear viscosity) to 37.2 Pa-s via changes 
in the suspension concentration, thus capturing low viscosities that have been historically very challenging to measure. These 
changes equate to an increase in the relative extensional viscosity of nearly 12,200 times at a volume fraction of just 0.027. We 
also find that interactions between the wires and the necessary polymer additive affect the rheology strongly. Polymer-induced 
elasticity shows a reduction as the buffer relaxation time falls from 819 to 59 ps above a critical particle concentration. The results 
and technique presented here should aid in the future formulation of these promising nanowire suspensions and their efficient 
application as inks and coatings. 


1 Introduction 

Copper nanowires are novel building blocks in nanotechnol¬ 
ogy, exhibiting a wide range applications in electronics, field 
emission transistors, heat management, and electrochemical 
detectors^. Networks of thin and long copper nanowires 
(CuNWs) (see Fig. |T|) demonstrate very good electron trans¬ 
port in sparse concentrations after deposition and low temper¬ 
ature annealing (100 to 200°C — permitting the use of flexible 
plastic substrates). This leaves large spaces between the wires 
for the transmittance of light, spanning UV to infra-red wave¬ 
lengths. Simultaneously, they also display high flexibly and 
stability-. These attributes make them appealing in the pro¬ 
duction of fully transparent and flexible conductive electrodes. 
Thus, they have been exploited in the development of devices 
such as light-emitting diodes, photovoltaic cells, stretchable 
conductors, photodetectors, solar control windows, and touch 
screens, to name a few— 

The utilisation of these particles is most likely to take the 
form of high-precision deposition or thin uniform coatings via 
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a deposited suspension. Thus, there is a need for inks and 
coatings that are suitable for inkjetting, spraying, and drawing 
(among others) to be formulated, requiring a balance between 
the bulk fluid properties, method of delivery, and the target 
substrate. These application processes are dominated by ex¬ 
tensional flow, complicating the development of suitable fluids 
for three reasons. First, these extensional processes predomi¬ 
nantly require low-viscosity fluids to function-’-’— (the buffer 
used in this study was only 1 mPa-s). Second, suspensions of 
CuNWs are likely to be highly non-NewtonianJi— ’— . Third, 
these particles require polymer additives in order to stay sus¬ 
pended, lest they aggregate and settle. Hence, phenomena 
observed in polymer solutions in extensional flow such as 
long lived filaments and beads-on-a-string— must also be ad¬ 
dressed. These challenges are compounded by the dearth of 
experimental extensional-flow methods that can handle low- 
viscosity fluids and quantify these surface-tension driven sys¬ 
tems. 

Inkjet technology generates extensional flows using low- 
viscosity fluids and it has been proposed that it can be used 
to characterise them—. However, inkjet systems have their 
limitations. Fromm (1984)— states that the critical properties 
that determine whether or not a particular fluid is suitable for 
inkjetting are the viscosity, density, surface tension, and the 
nozzle diameter. Given typical nozzle diameters of 10 and 
40 pm, ink viscosities of 2-20 mPa-s and 5-50 mPa-s are rec¬ 
ommended, respectively, as evidenced by recent particle sus¬ 
pension ink work—’—’—’—’—’—’—. These viscosity ranges 
are not only narrow but also range below the lower limit of 
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Fig. 1 Network-forming nature of long and thin CuNWs 
(approximately 27.6 nm in diameter, 62.36 pm long, and possessing 
an aspect ratio of 2200 on average). A vast number of connections 
that is a result of annealing can be seen while a large area for light to 
penetrate is maintained. 


technologies like CaBER, an important consideration for char¬ 
acterising the complex low-viscosity buffers and dilute sus¬ 
pensions in this study. This point is further complicated here 
since the viscosity of these high-aspect ratio particle suspen¬ 
sions is anticipated to be highly sensitive to volume fraction, 
a potential hurdle in ink formulation when balancing conduc¬ 
tivity against light transmittance. Further, it is recommended 
that the nozzle diameter is at least 100 times greater than the 
suspended particles^, limiting the range of sizes that can be 
handled by any one system - an issue likely complicated by 
these very high-aspect ratio particles, by which logic would 
potentially require nozzles of 6236 pm. Further compounding 
the use of inkjets in the analysis of these long and thin wires 
are the complex flows and large shear rates generated within 
the nozzle itself that can be 2 x 10 4 to 2 x 10 6 s" 1 — , which 
may cause damage to NWs as they do to delicate proteins^’—. 
Moreover, if thermal inkjet technology (as opposed to piezo) 
were used then the raised temperatures of around 300°C ex¬ 
ceed the annealing temperature for these particles, probably 
damaging them and possibly causing nozzle clogging. Ex¬ 
acerbating all this is the fact that every inkjet system has its 
own characteristics that produce system-specific results that 
can vary considerably. Though inkjet technology may be use¬ 
ful for delivering these suspensions in certain formulations, it 
is unlikely to be suitable for characterising them or complex 
buffers. Hence, finding a workable recipe to make inks from 
these NWs would otherwise be trial and error. 

Although accepted shear rheometry techniques have long 


been available 2 -, established extensional measurement meth¬ 
ods are comparatively recent—’—. Capillary-breakup exten¬ 
sional rheometry (CaBER) involves rapidly extending a drop 
of sample fluid between two end-plates to a finite distance, 
creating a liquid filament with a critical aspect ratio that then 
thins due to surface tension. The thinning evolution is de¬ 
scribed by the balance of surface tension against inertial and 
viscous stresses caused by the extensional flow. Through this 
balance it is possible to extract the extensional viscosity by 
measuring the decay of the neck radius as a function of tinted. 
This technique has been used for highly viscous samples^ but 
encounters problems with low-viscosity fluids. Rapid stop¬ 
page of the end plates causes vibrations to propagate through 
the fluid filament, leading to inertial instabilities which are not 
accounted for in the stress balance. These instabilities dissi¬ 
pate and can be ignored in higher viscosity samples or in sam¬ 
ples with elasticity, where long-lived cylindrical polymeric fil¬ 
aments allow the use of a simple stress balance to extract the 
viscoelastic stresses^. 

Recent advancements of the CaBER concept provide fur¬ 
ther access to low-viscosity fluids. The Cambridge Trimaster 
increases the end-plate retraction speed to overcome inertial 
instabilities and allow Newtonian and weakly elastic fluids of 
10 mPa-s to be tested—’—. Another is the Slow-Retraction- 
Method—, in which a liquid filament is extended between 
two end-plates to a stable aspect-ratio, followed by a slow 
extension that pulls the filament beyond the critical aspect- 
ratio to initiate break-up. Electrowetting forces can also be 
used to pull a fluid filament to initiate capillary thinning—. 
A liquid to liquid annular flow microfluidic device has also 
been developed—. The recently developed Rayleigh Ohne- 
sorge Jetting Extensional Rheometry— technique involves 
stroboscopic imaging of a continuously flowing jet with con¬ 
trolled surface perturbations and has allowed access to very 
thin weakly viscoelastic fluids, though higher sample volumes 
would make such a device unsuitable for this study. 

We use a novel variant of the CaBER technique that uses 
surface acoustic waves (SAW) to create and stabilise liquid 
bridges of very low viscosities; distinct advantages of this 
technique make it suitable for the development and analysis 
of these NW suspensions as inks and coatings. Firstly, it can 
process Newtonian fluids possessing shear viscosities from 27 
Pa-s down to 0.96 mPa-s, giving it the range to encompass both 
viscous and characteristically difficult low-viscosity fluids^. 
Moreover, despite delivering large bulk mechanical energy the 
shear amplitudes created by high frequency SAW are so small 
that they can manipulate delicate particles, such as DNA, pro¬ 
teins and cells, without causing damage^’—’—’—. Further, 
SAWs have also been employed to pattern liquid crystals^ 
and nanotubes^ without damage. Additionally, heating ef¬ 
fects for the short bursts required to jet a droplet of water or 
glycerol have been determined to be less than 0.02 and 1° C, 


2 | Journal Name, 2010, [volfl -[7] 


This journal is ©The Royal Society of Chemistry [year] 










10 - 


8x1 O' 4 


6x1 O' 4 


4x1 O' 4 


2x1 O' 4 


b 


1 1 1 1 1 1 

[ buffer 

-iiiiiiiiiiiiiiiiii 

m 7.5x10" 5 

- 

; 8.3x10 -5 


9.2x10" 5 

- All 8x10 -4 

°A-1.66x10" 4 


1.67x10" 3 2.5x10 -3 

: //" *\ ■ D o m 3 \ 10 ’ 3 4.15x1 O' 3 ' 

; 4.15X10- 4 \ 3x10 -4 1-3 l a 1 ; 

- 2.76x1 0~ H 

_i_i_i_i_i_i_ 

3x10"° ° 


10- 3 


2x1 O' 2 


3x1 O' 2 


4x1 O' 2 


5x1 O' 3 


Fig. 2 (a) Filament break-up data clearly showing the reduction in polymeric effect with an increase in the CuNWs:polymer ratio. The red 
and green lines indicate the polymeric regions used calculate the relaxation times A e for the buffer and 4> = 4.15 x 1CT 3 , respectively. Break-up 
above R* = 1/2 shows Newtonian-like behaviour, allowing extensional viscosity extraction as described in Section 2.1. The optical lower 
limit in these experiments occurs at 4 pm. (b) Relaxation times showing a stark drop from 819 to 248 ps at a CuNW volume fraction of 4> = 
1.18 x 1CT 4 , indicating a drop in the amount of free polymer in the solution. A more gently approached reduction to 59 ps occurs at 
4.15 x 1CT 3 . Grey annotations denote volume fractions. 


respectively — favourable for rheological analysis and well 
below the annealing temperatures of these metallic particles. 
Furthermore, there is no nozzle that can be clogged or that can 
restrict particle size. Finally, the fact that the technique only 
requires micro-droplets lends itself to the analysis of novel flu¬ 
ids such as biological samples that are difficult to produce in 
large volumes^ (another microfluidic rheometric technique 
has also been recently proposed to handle small sample vol¬ 
umes^). This is an important advantage since it is currently 
difficult to synthesise large quantities of CuNW particles. 

Hence, we report changes in the extensional viscosity 
as a function of the particle concentration from dilute to 
highly concentrated and viscous suspensions. Additionally, 
the extensional properties of the polymer-doped aqueous 
buffer were also analysed, with unexpected phenomenological 
changes as the CuNW:polymer ratio was changed. The avail¬ 
ability of this data and technique is anticipated to assist in the 
development and application of these promising particles. 

2 Results and discussion 

2.1 Polymeric buffer analysis 

Phase separation and settling due to depletion interaction has 
been reported in previous experiments with suspensions of 
particles in polymeric buffer—. However, we observed no set¬ 
tling in our samples. This suggests that in our case polymers 


may absorb onto wires and serve to repel other coated wires, 
preventing aggregation and stabilising these suspensions. 

Capillary thinning of pure polymer solutions has been 
widely studied— It is known that when a liquid 
bridge starts thinning, polymeric stresses are negligible, un¬ 
til polymer molecules have stretched sufficiently in the exten¬ 
sional flow at the necking plane. Molecular stretching leads to 
strain-hardening of the polymeric stresses and at some point, 
these stresses begin to dominate over the viscous stresses from 
the rest of the fluid. At that point, the decrease in filament ra¬ 
dius changes from an approximately power-law like decay in 
low-viscosity solvents to an exponential decay. The relaxation 
time of this exponential decay is used as an indicator of the 
elasticity induced by the polymeric component. 

We are interested in extracting the viscosity of the parti¬ 
cle suspension alone without the interference from the poly¬ 
meric component. A reduction in polymeric effects, quanti¬ 
fied by the relaxation time A e , is seen with an increase in the 
CuNW:polymer ratio (see Fig. 0. The relaxation time drops 
from 819 ps for pure buffer to 59 ps at a volume fraction of 
4.15 x 1CT 3 . We reason that as polymers bond to the wires 
they lose their ability to stretch under flow and the addition of 
more wires leaves fewer polymer molecules free to influence 
bulk fluid behaviour. We then use the CuNW:polymer ratio 
at which A e is the smallest, while still maintaining the sus¬ 
pension, to minimise interference from the polymer, in either 
shear or extensional viscosity measurements. 
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Fig. 3 (a) Relative extensional f\ re i and shear viscosity r\ re x 
showing similar behavioural trends as a function of the volume 
fraction 0. The shear viscosity values correspond to comparable 
extensional and shear strain rates. It should be noted that Sec l3.2l 
describes how the extensional viscosities seen here were extracted, 
(b) Relative shear viscosity against volume fraction for shear rates 7 
of 20, 40, 80, 160, and 200 s' 1 , (c) Relative shear viscosity against 
shear rate for different volume fractions, displaying pronounced 
shear thinning. 


2.2 Viscosity measurements 

Fig. [3^ shows the progressive change in the relative exten¬ 
sional viscosity f\ re \ - fj/3rj s with nanowire volume fraction, 
where r] s =1.01 mPa-s is the shear viscosity of the buffer with 
the polymeric solute. A comparable growth in shear viscosi¬ 
ties is also observed. It must be noted that each data point 
in Fig. [3^ is at a different strain-rate. In capillary break-up 
experiments, the strain-rate cannot be set a priori but is self- 
determined by the fluid. In general, the average extensional 
strain-rate experienced by a sample prior to viscosity extrac¬ 
tion decreases as volume fraction increases and viscosity in¬ 
creases. In shear rheometers on the other hand the strain rate 
can be pre-set. The shear data in Fig. [3 \l were obtained at 
strain-rates comparable to the corresponding extensional flow 
experiment at the same volume fraction. The shear viscosity 
data obtained at fixed strain-rates are shown in Figs. 0> and 


c. Since the polymeric component is expected to shear thin, 
and at any given volume fraction, the zero-shear-rate viscosity 
contribution of the polymer is small compared to the nanowire 
contribution (as demonstrated by the lack of significant expo¬ 
nential decay at small times during capillary break-up), the 
shear data in Figs. 0) and c can also be considered to be free 
of polymeric interference. 


Acute viscosity sensitivity to volume fraction is expected 
for these CuNW suspensions since these wires have the poten¬ 
tial to directly interact with large numbers of particles along 
their lengths, and thus quick growth in the effects of inter¬ 
particle interaction relative to increased volume fraction was 
anticipated. In fact, the extensional viscosity more than dou¬ 
bled the buffer viscosity at a volume fraction of 1 x IO -3 . At 
4> = 1.3 x IO -2 , an extensional viscosity of 265 mPa-s was 
reached, well in excess of the typical inkjet system capabil¬ 
ity discussed in Section 1. Further, the relative extensional 
viscosity nearly rose to 12,200 times (77 = 37.2 Pa-s) at only 
0 = 2.7 x 1CT 2 . The steep gradient at this point indicates that 
the maximum particle packing fraction was being approached 
even at this low volume fraction, which may be expected for 
such high aspect ratio particles^. 


It can be seen in Fig. [3 \l that the shear viscosity overtakes 
the extensional viscosity in magnitude at about <j> = 1.7 x 1CT 3 . 
To our best knowledge, this has not been observed before. 
The reasons for this effect are not clear at present and may 
be attributed to the difference in the interplay between flow 
and particle orientation. The bare nanowires have an aspect- 
ratio of approximately A = 2 x 10 3 . Even after accounting 
for an absorbed polymer layer, the aspect-ratio A * O(10 3 ). 
The volume fraction for the isotropic-nematic transition is ex¬ 
pected to be around 1/A » IO -3 . Most of the samples in 
Fig. 0 are well above this concentration. Extensional rhe¬ 
ology of liquid-crystalline materials is relatively unexplored 
and the few studies thus far have focused on high-temperature 
liquid-crystalline polymer melts^’— . Extensional viscosi¬ 
ties larger than shear viscosities have been observed in liquid- 
crystalline carbon-nanofibre suspensions, but at a single vol¬ 
ume fraction—. CuNW suspensions could serve as model 
room-temperature liquid crystalline materials for more de¬ 
tailed explorations for understanding the effect of particle con¬ 
centration on behaviour in extensional and shear flows. On the 
applied front, the data in Fig. [3 \l show that to a large extent, 
shear viscosity can provide a good estimate of extensional vis¬ 
cosity for use in designing inkjet applications of CuNW sus¬ 
pensions. 
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3 Experiments 

3.1 Particle and buffer preparation 

CuNWs 25-30 nm in diameter and 50-65 pm long were syn¬ 
thesised via a hydrothermal method—. 900 mg hexadecy- 
lamine (HDA)(Sigma-Aldrich) was dissolved in 50 ml water 
in a closed 100 ml screw cap borosilicate glass bottle, and 
magnetically stirred at 1000 rpm at 100°C for 1 hour in an 
oil bath. To this, 100 mg CuC 1 2 2H 2 0 was added and stirred 
for another 30 minutes, after which 500 mg of a or /3 glucose 
(Merck) was added and stirred at 400 rpm for 6 hours. The 
reaction solution takes on the appearance of a creamy white 
liquid with just the HDA after 1 hour. After the addition of 
CuC 1 2 2H 2 0 the solution turns sky blue after 30 minutes. Fi¬ 
nally, after addition of glucose, the solution starts turning pale 
brown at first, then gradually dark brown, then reddish brown. 
The reaction is stopped after 6 hours and removed from the oil 
bath allowing it to cool for 10 minutes, after which it is trans¬ 
ferred into 50 ml vials and centrifuged at 6500 rpm for 5 min¬ 
utes. The CuNWs are collected at the bottom of the centrifuge 
tube, which can be easily separated by carefully decanting the 
waste supernatant. The pellet is carefully rinsed a few times 
with Milli-Q water at room temperature to remove any loosely 
bound HDA remaining. 

Hydroxypropyl cellulose (HPC)(Sigma Aldrich, 20 mesh, 
Mw = 1,000,000) was the polymer used as a dispersant to help 
suspend the CuNWs and prevent aggregation. A 2% wt so¬ 
lution of HPC in ethanol was made from of which a 0.6 ml 
aliquot was sperated into a test tube and diluted with Milli-Q 
water to make a 2.2 ml solution. This aqueous solution was 
used as the buffer, and the same amount of HPC solution was 
used with increasing volume fractions of CuNWs to find the 
optimum CuNW:polymer ratio through extensional analysis. 

3.2 Acoustically-driven microfluidic extensional rheome- 
try 

We overcame the challenges discussed in Section 1 by using 
an experimental technique developed by Bhattacharjee et al ., 
explained in detail elsewhere^’—. Briefly, we used a 20 MHz 
waveform generator (33220A, Agilent) to trigger a second sig¬ 
nal generator (WF1966, NF Corporation), which produces a 
1.5 ms burst. The second generator outputs a sinusoidal signal 
to an RF power amplifier (411 LA, ENI), resulting in a fixed 
amplitude and frequency signal of approximately 36.7 MHz 
that is then supplied to a piezoelectric SAW chip. This de¬ 
vice comprises of a pair of curved inter-digitated transducers 
bonded to a piezoelectric lithium niobate substrate that focus 
Rayleigh waves towards a central target where a sessile droplet 
in placed (1 pi; approximately 1 mm dia - Fig 01) The SAW 
burst, focused beneath the droplet, concentrates considerable 
inertial energy within it, thus forcing a jet to form—. The SAW 



'1 mm 


Fig. 4 A montage showing the SAW-actuated jetting, resulting 
from a 1.5 ms SAW burst, and subsequent surface tension-driven 
thinning of pure water at a temperature of 22°C (0.96 mPa-s). 


energy then stabilises the jet as it impacts onto an opposing 
surface and establishes a liquid bridge of a critical aspect ratio 
that will result in break-up. Upon relaxation of the SAW burst, 
the filament is allowed to thin under surface tension, the decay 
of which is recorded using a high speed camera (Photron S A5; 
62000 fps; image size: 1.35 mmx2.14 mm (192x304 pixels) 
with a long-distance video microscope attachment (K2/SC, In¬ 
finity). The event is lit by a single LED lamp placed beyond 
the fluid filament. The radius of the neck in each image frame 
is extracted using standard image-analysis techniques. 

The method for extracting extensional viscosity for this 
technique is detailed elsewhere^. In short, the aforemen¬ 
tioned difficulties of extracting data from thin cylindrical 
threads are addressed by beginning the measurement of fila¬ 
ment thinning at a larger filament diameter, but not so large 
that initial transients affect data — 50 pixels (0.352 ± 0.007 
mm). A series of filament break-up data was acquired for a 
range of Newtonian fluids that represent viscosity standards, 
from low-viscosity to viscous samples. Break-up time and 
radius were re-scaled by the Rayleigh time-scale t * = £/tr 
and the initial radius R* = R/Ro, respectively, where tr = 
\[pR o/ a, pis the density, Rq is the initial radius, and a is the 
surface tension. When experimental conditions such as sam¬ 
ple size were held constant, it was determined that the time 
t\j 2 for R* to decrease to a value of half was primarily a func¬ 
tion of the Ohnesorge number Oh= fj/3\/paRo , where ij is the 
extensional viscosity, which is 3 times the shear viscosity p for 
Newtonian fluids. A calibration curve through t\j 2 - versus - 
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Oh data for the Newtonian viscosity standards was obtained. 
With this calibration curve, the Oh and hence fj of a test fluid 
sample with negligible elastic effects can be extracted from 
the observed t\j 2 measured. This method for extracting ex- 
tensional viscosity is valid as long as viscoelastic effects due 
to the dissolved polymer are not significant, since the origi¬ 
nal calibration of dimensionless half-time against Ohnesorge 
number was only performed with Newtonian viscosity stan¬ 
dards. 

3.3 Characterisation of polymeric filaments 

The addition of polymer to the suspension buffer necessi¬ 
tates quantification of viscoelastic properties. Polymer chains 
stretch under extensional flow, retarding the rate at which a 
fluid filament breaks up. These elastic stresses can domi¬ 
nate towards the later stages of break-up, resulting in long- 
lived filaments and can be characterised by a relaxation time 
A e = -f/31n( J R(t)/i?o)“ , “ • Although in reality there may 
actually be more than one relaxation mode, it has been found 
that the slowest mode dominates in extensional flow, thus a 
single relaxation time is commonly used to characterise vis¬ 
coelastic stress experimentally and theoretically—. 

3.4 Particle imaging 

High resolution images of the CuNWs were captured using 
field-effect scanning electron microscopy (Nova NanoSEM 
450, FEI, Oregon, USA). 

3.5 Shear rheometry 

A cone (60 mm/l°, gap distance 0.041 mm) and plate rheome¬ 
ter was used to perform the shear viscosity analysis (HAAKE 
MARS III, Thermo Fisher Scientific, Massachusetts, USA). 

4 Conclusions 

The capabilities of SAW-driven microfluidic extensional 
rheometery enabled the extensional analysis of NW suspen¬ 
sions, overcoming the challenges that low-viscosities, high- 
aspect ratio particles, and small sample volumes present. 

The reduction of polymeric-induced elasticity with the ad¬ 
dition of CuNWs, indicated by the relaxation time decreas¬ 
ing from 819 down to 59 ps, suggests that polymer molecules 
bond to the CuNWs. Beyond a certain CuNW:polymer ra¬ 
tio, few free polymer molecules remain to influence bulk 
properties. Under these conditions, filament thinning was 
initially not significantly affected by the dissolved polymer, 
and Newtonian-like behaviour was observed at small strains. 
Hence, utilising this optimal ratio enabled the creation of 


sedimentation-free suspensions that did not exhibit compli¬ 
cating polymeric behaviours, thus allowing the observation 
of changes in viscosity that arise primarily as a function of 
CuNW volume fraction. 

It was demonstrated that suspension viscosity of these long 
and thin wires are highly sensitive to volume fraction, with 
increases of 4 orders of magnitude in their relative extensional 
and shear viscosities, at a volume fraction of just 2.7 x 1CT 2 . 

In summary, it is anticipated that the reported technique and 
viscosity data, particularly that of problematic low-viscosity 
fluids in extensional flow, will assist in the formulation of 
workable inks and coatings where controlling polymeric ex¬ 
pression and extensional viscosity is key, thus expediting the 
application of these promising particles for the next generation 
of electronics. 
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Fig. 5 For Table of Contents Only. Novel experiments enabled the 
extraction of lacking extensional data on promising copper nanowire 
suspensions and their complex buffer, where low viscosities and 
long ID particles pose significant challenges for alternative 
technologies. 
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